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Abstract 

This paper reports on the experimental implementation of the quantum baker's map via a three 
bit nuclear magnetic resonance (NMR) quantum information processor. The experiments tested 
the sensitivity of the quantum chaotic map to perturbations. In the first experiment, the map was 
iterated forward and then backwards to provide benchmarks for intrinsic errors and decoherence. 
In the second set of experiments, the least significant qubit was perturbed in between the iterations 
to test the sensitivity of the quantum chaotic map to applied perturbations. These experiments 
are used to investigate previous predicted properties of quantum chaotic dynamics. 

PACS numbers 03.67.-a, 03.67.Lx, 05.45.-a 



Chaos is a phenomenon in which nonlinear dy- 
namical systems exhibit heightened sensitivity to 
small perturbations H1H§||. 
The study of chaos is computationally intensive. 
When direct experiments are not available, com- 
puters can be used to simulate chaotic dynamics 
and to calculate the effect of perturbations. But 
when the chaotic system is quantum-mechanical, 
simulating its dynamics on a classical computer 
is notoriously difficult: the computational com- 
plexity of the calculation rises exponentially with 
the number of degrees of freedom of the simu- 



lated quantum chaotic system and with the ac- 
curacy to which the simulation is to take place 
[f| || [10]. If one can simulate quantum chaos on 
a quantum computer, by contrast, the computa- 
tional complexity rises only small polyno- 
mial in the number of degrees of freedom and in 
the accuracy [10] [11] [12|. Consequently, quan- 
tum computation represents a potentially pow- 
erful technique for investigating quantum chaos. 
This paper reports on an experimental demon- 
stration of a quantum chaotic map. 

Small perturbations to the initial state of a 
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classical chaotic system typically lead to large 
changes in behavior. Two states that are ini- 
tially close are driven apart at a rate governed by 
the positive Lyapunov exponents of the chaotic 
dynamics. By contrast, quantum dynamics, 
whether regular or chaotic, preserves the over- 
lap between quantum states and does not drive 
them apart. Nonetheless, quantum chaos can be 
characterized by the sensitivity of the time evo- 
lution of states to small changes in the Hamilto- 
nian that governs the chaotic dynamics || [|J §1] 

III 

Peres f|j|7] noted that under quantum chaotic 
dynamics with a slight perturbation, a state 
moves apart at an exponential rate from the 
same state evolving under the unperturbed dy- 
namics. Schack and Caves [|| characterized 
chaotic dynamics (both classical and quantum) 
in terms of the exponential growth of the infor- 
mation required to specify a state that evolves 
according to a perturbed version of the dynam- 
ics, a phenomenon they termed 'hypersensitivity 
to perturbation.' Zurek and Paz [13] conjectured 
that quantum chaotic systems are sensitive to 
weak interactions with the environment, which 
cause such systems to produce information at 
a rate equal to their Kolmogorov-Sinai entropy. 



Ballentine and Zibin [14] compare quantum and 
classical chaotic systems by measuring the accu- 
racy with which the system returns to its initial 
state under time reversal after a perturbation 
of varying strength. The implementation of the 
baker's map reported here allows for the experi- 
mental investigation of quantum chaotic dynam- 
ics, specifically its sensitivity to perturbations as 
described by the above models. 

The quantum information processor used to 
simulate perturbed quantum chaotic maps is an 
NMR quantum information processor |16| Jl7|]. 
The number of quantum bits used (three) is suf- 



ficiently small that the precision of the quantum 
computation can be checked on a classical com- 
puter. Of course, the small number of qubits 
means that the simulation could have been per- 
formed on a classical computer. The goal of 
the research reported here was to actually sim- 
ulate quantum chaos on a quantum information 
processor (as opposed to using a classical com- 
puter to simulate a quantum information pro- 
cessor simulating quantum chaos |ll|). We also 
explore techniques and problems for scaling up 
such simulations to allow for a thorough inves- 
tigation of quantum chaotic dynamics of high- 
dimensional quantum systems. A further goal of 
experimental investigations is to identify practi- 
cal experimental signatures of quantum chaos. 

The quantum chaotic map investigated is the 
quantum baker's map [|l^] [19]. The classical 
baker's map acts on the unit square in phase 
space as follows: The baker's map first stretches 
phase space to twice its length, while squeezing 
it to half its height. Then, the map cuts phase 
space in half vertically and stacks the right por- 
tion on top of the left portion, similar to the way 
a baker kneads dough. Because of the stretching 
and the cut, the baker's map is fully chaotic and 
has two Lyapunov exponents, ±ln2. Balazs and 



Voros [18 1 presented a quantized version of the 
baker's map that reproduces the behavior of the 
classical map in the limit 7l — > 0. The quantum 
baker's map (QB) is a simple unitary operator 
which consists of a quantum Fourier transform 
(QFT) on half of the Hilbert space followed by 
an inverse QFT on the whole Hilbert space. 

The quantum baker's map can be expressed as 
a sequence of two basic unitary operations, the 
Hadamard gate, Hj, operating on spin j and the 
conditional phase gate, Bjk, operating on spins j 



and k. Schack [12] utilized this decomposition to 
develop an algorithm for simulating the quantum 
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baker's map on a quantum computer. The three 
qubit version of the quantum baker's map is 

Swapi 3 H 3 Bl 3 Bl 3 H 2 Bl 2 H 1 Swap23H 3 B23H 2 

(1) 

where Swapj^ is a swap gate between bits j and 
k. This gate sequence is a QFT on qubits 2 and 
3 [20 1, followed by an inverse QFT on all three 



qubits. 

We note that according to a well known con- 
jecture ||[21]Q, the energy eigenvalues of a 
quantum system with classically chaotic dynam- 
ics are expected to have the statistical proper- 
ties of a Gaussian random matrix. A quantum 
baker's map whose Hilbert space dimension is a 
power of two (which is the case for our implemen- 
tation) does not exhibit this spectral distribution 
Nevertheless, it does show interesting be- 
havior due to perturbations. 

The approach taken in this implementation 
of the quantum baker's map is to apply a 
perturbation between forward and backward 
map iterations. In general, the final state 
with perturbation can be expressed as = 

E k (QB^A k (QB)p inmal (QB)^Al(QB), whereas 
the final state without perturbation is pj = 
(QB) J <(QB)p initia i (QB)^(QB) In the absence of 
experimental errors Pf = pinitial. 

An interesting non-unitary perturbation is the 
application of a magnetic field gradient dephas- 
ing the least significant (third) bit. The Hamilto- 



nian of a gradient on spin j is H 3 g (z) = e >~dT z ~ ; 
and it acts clS el rotation of varying magnitude 
across the sample. This perturbation looks like 
decoherence when measuring the signal from the 
entire sample, effectively tracing over position, 
z. 

An interesting type of unitary perturbation 
is one with a simple definition in phase space, 



such as a displacement in position. The quan- 
tum baker's map, as defined by equation (2), 
associates the discretized position basis on the 



unit square [18] with the computational basis. 
The smallest possible position perturbation in 
a Hilbert space of dimension N corresponds to 
the addition or subtraction of l(modN) which, 
for example, would take the state [000) to [001). 
This perturbation can be implemented on n 
qubits by applying a (j — l)x controlled-NOT 
to bit j, where the controls are all bits less sig- 
nificant then j, as j is indexed from n to 1. Us- 



ing geometric algebra 1 22 ] this sequence for three 
bits can be written as follows: 



5 *f(i-«i-«2-«»2 



(2) 



Simulated data in the inset of figure (1) shows 
that, as expected, the baker's map is very sensi- 
tive to such a perturbation. 

A variation of this perturbation with a simple 
interpretation in Hilbert space is an x rotation 
on one bit. The sensitivity of the baker's map 
to the a x perturbation may be carefully exam- 
ined experimentally by controlling the angle of 
the rotation. In our experiments we study only 
one initial state, the [000} psuedo-pure state. In 
general, the perturbed dynamics of a single state 
may not be sufficient to completely characterize 
the behavior of a map. Here, we take advantage 
of the fact that the overlap Trace{pfp'^) for the 
[000) initial state approaches the average of over- 
laps for a complete set of orthogonal initial states 
in the limit of small perturbation, as shown in 
figure (1). Thus the overlap, Trace(pfp'j) for the 
initial state [000) approaches the fidelity p3] [24] 



of the operation: baker's map, rotation pertur- 
bation, inverse baker's map. 

The three bit quantum baker's map was im- 
plemented via NMR using the three carbon-13 
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spins of an alanine sample. The resonant fre- 
quency of carbon-13 on a 300MHZ spectrometer 
is approximately 75.468MHz. Frequency differ- 
ences between the spins are 9456.5Hz between 
spins 1 and 2, 2594.3Hz between spins 2 and 
3, and 12050.8Hz between spins 1 and 3. Cou- 
pling constants between the three spins are J\2 
= 54.2Hz, J 23 = 35.1Hz, and J 13 = -1.2Hz. Re- 
laxation time T\ for the three carbon spins in 
alanine are all longer than 1.5s while the T 2 re- 
laxation times are longer than 400ms. The pulse 
sequences for realizing the Hj and Bj^ gates and 
the implementation of the QFT as a sequence 



of these gates, can be found in [25]. The pulse 
sequence for the complete quantum baker's map 
was compressed by relabeling bits instead of per- 
forming the swaps explicitly. Readout was done 
using quantum state tomography as described in 



Figure 1: Sensitivity of baker's map to unitary 
perturbations. Solid line shows the theoretical 
overlap Trace(pfp'j) versus the angle of the per- 
turbation rotation on the third bit for the initial 
state |000). The dashed curve shows the final 
state overlap for the same perturbation averaged 
over a complete set of orthonormal initial states. 
For small perturbations the results for the initial 
state 1 000) are very similar to that of a complete 
set of initial states. Stars represent the mea- 
sured overlap Trace(gfg'j) for the initial state 
1 000). The inset plot is the overlap versus the 
size of the position perturbation (add 1 through 
add 4). The baker's map is much more sensi- 
tive to a position space perturbation then to a 
rotation perturbation on the least significant bit. 



To measure the accuracy with which the trans- 
formations were performed we used the correla- 
tion measure introduced in |24|] which is appro- 
priate for almost fully mixed density matrices: 



C 



Tr ( Qtheory Qexp) 



Tr(g? 



expJ 



eory> 



Tr{g 2 



exp/ 



T 7 '(S 'initial) 



(3) 

Here, g is the measured deviation density ma- 
trix. This does not include the large identity 
term in the actual high temperature liquid state, 
but does include a small amount of the identity 
operator necessary for reconstructing a positive 
state operator. The amount of identity included 
is fixed for each set of experiments. If the the- 
oretical and experimental deviation density ma- 
trices were correlated C = 1, if they were uncor- 
rected C = and if they were anti-correlated 
C = -l. 

The first experiment consisted of two itera- 
tions of the quantum baker's map, a forward it- 
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eration followed by the inverse of the map, start- 
ing from the |000) pseudo-pure state. The at- 
tenuated correlation, C, of the implementation 
of the forward map is .76, and for the forward 
followed by the inverse, .56. The first term in 
the attenuated correlation measures the correla- 
tion between theoretical and experimental den- 
sity matrices without accounting for reduction 
in signal over the course of the experiment. This 
gives a crude measure the maps' accuracy in the 
absence of decoherence. For the forward map 
this unattenuated correlation is .93 and for the 
forward followed by the inverse, .90. Since the 
experiment was done on the |000) pseudo-pure 
state, we expect the final state of the system to 
be that same state. The density matrix of the 
spin system after the forward and inverse itera- 
tion is shown in figure (2). 

Another set of experiments were done to ex- 
plore the dynamics of the baker's map under per- 
turbations as described above. The first pertur- 
bation experimentally tested consisted of a de- 
phasing gradient on the third, least significant, 
bit. We measured C = .78 for the state after the 
gradient and C = .65 for the state after the in- 
verse map. For these correlation values the loss 
of magnetization due to the gradient is taken into 
account in the normalization of Qtnitial- 

The fi- 
nal state (after the gradient and inverse map) 
shows the generation of one bit of entropy as 
seen by the equilibration of the |000) and |001) 
populations, as displayed in figure (2). This is 
consistent with the Paz-Zurek model for the ef- 
fect of decoherence on a quantum chaotic map: 
when decohered, the map produces one bit of 
entropy per iteration, an amount equal to the 
Kolmogorov-Sinai entropy of the map. 

In a final set of experiments rotational per- 
turbations of 7r/32, 7r/16, 7r/8, and 7r/4 on the 
third (least significant) bit were applied between 




Figure 2: Real part of experimental density ma- 
trices. The top left hand corner shows the in- 
put 1 000) psuedo-pure state. One iteration of 
the baker's map on this state leads to the center 
density matrix. Then the inverse map is applied 
(top right figure) bringing the bits back to the 
input |000) state. In the second experiment, a 
gradient is applied to the third, least significant, 
bit after the forward iteration of the map (bot- 
tom left figure) followed by the application of the 
inverse map (bottom right figure). 
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the forward and inverse iterations of the baker's 
map. For these experiments C ranged between 
.52 and .53. In figure (1) we plot the overlap, 
Trace(gfQ'f), of the experimental perturbed and 
unperturbed density matrices and compare it to 
theoretical predictions. 

In conclusion, we describe the implementation 
of a chaotic map on a quantum system, a three 
qubit quantum information processor. In addi- 
tion, we have explored two perturbations and ex- 
amined their effects on the dynamics of the map. 
Experiments such as these establish a foundation 
for further experimental investigations of quan- 
tum chaotic dynamics and the exploration of 
suggested theoretical approaches. For example, 
the hypersensitivity to perturbation suggested 
by Schack and Caves, should be evident even 
on a small Hilbert space, with only a few iter- 
ations of a chaotic map [11|. To eventually ob- 



serve a characteristic such as the Peres criterion 
H will require a much larger Hilbert space and 
many more iterations of the map. We believe the 
experiments performed are a first step towards 
a more thorough experimental investigation of 
these questions. 

The authors thank M.A. Pravia and E.M. 
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